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Bose-Einstein Condensates 



OO 
0> 



Elena V. Goldstein and Pierre Meystre 
Optical Sciences Center, University of Arizona, Tucson, AZ 85721 

We present an exact quantum mechanical analysis of collinear four-wave mixing in a multicompo- 
nent Bose-Einstein condensate consisting of sodium atoms in the F = 1 ground state. Technically, 
this is achieved by taking advantage of the conservation laws of the system to represent its Hamil- 
tonian in terms of angular momentum operators. We discuss explicitly the build-up of matter-wave 
side-modes from noise, as well as the correlations between these modes. We show the appearance of 
a strong quantum entanglement between hyperfine states. We also demonstrate that for finite atomic 
numbers, the system exhibits periodic collapses and revivals in the exchange of atoms between dif- 
ferent spin states. 

PACS numbers: 03.75.-b 03.75.Fi 05.30.Jp 42.65.Hw 
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I. INTRODUCTION 

Low density Bose-Einstein condensates of alkali atoms 
are described to an excellent degree of approximation by 
a Gross-Pitaevskii nonlinear Schrodinger equation. Such 
equations are ubiquitous in many fields of physics, in- 
cluding nonlinear optics, hence it is not surprising that 
many of the concepts first developed in optics can read- 
ily be extended to condensates. A number of theoretical 
investigations along these lines have already been pre- 
sented, including the study of matter- wave solitons [Q, 
phase conjugation |2|,||], four-wave mixing pj, etc. The 
experimental verification of these predictions is also un- 
der way, with the first announcement of four- wave mixing 
in a rubidium condensate ||. 

A recent experimental development of considerable im- 
portance in this context is the demonstration of multi- 
component condensates, in particular of 23 Na conden- 
sates in far-off-resonance optical dipole traps ||. These 
systems permit to distinguish condensate modes in ways 
other than by their center-of-mass quantum numbers. 
This opens up new directions of research, such as the 
study of the stability and miscibility of quantum fluids, 
the analysis of pattern formation, the generation of ferro- 
magnetic and antiferromagnetic states, etc... [@,|]]4j|5|] . In 
the context of nonlinear atom optics, the coexistence of 
condensates with different magnetic quantum numbers is 
attractive in that it provides a way to perform four-wave 
mixing experiments in collinear geometries, with the con- 
siderable advantage of eliminating phase-matching limi- 
tations ||. 

It is now well understood that the matter-wave ana- 
log of nonlinear optical interactions is provided by colli- 
sions. In particular, two-body collisions in the shapeless 
approximation are mathematically equivalent to a local 
Kerr medium with instantaneous response. Hence, the 
kinds of wave-mixing phenomena that can take place in 
a condensate are largely dictated by the properties of 
these collisions. For example, recent experimental stud- 
ies of spin population dynamics in a 23 Na condensate 



have shown that as a result of spin changing collisions, 
a sample initially in the mp — state eventually winds 
up with almost equal populations of all hyperfine ground 
states The present paper exploits these collisions to 
generate four-wave mixing in 23 Na. It is organized as fol- 
lows: Section II describes our physical model, which con- 
siders degenerate backward phase conjugation in (small) 
condensates. Section III derives a representation of the 
problem in terms of angular momentum algebra. This 
leads to an exact diagonalization of the four-wave mixing 
Hamiltonian and the determination of all eigenstates and 
eigenenergies. These results are applied in Section IV to 
the analysis of the dynamics of exchange of population 
between different spin states. The quantum correlations 
between these states are also discussed. Finally, Section 
V is a summary and outlook. 

The phase conjugation of atomic waves has previ- 
ously been discussed in a situation where the matter- 
wave modes involve the same electronic state, but dif- 
ferent center-of-mass components, and using the unde- 
pleted pump approximation @. Mode separation was 
achieved via Bragg scattering off an optical field. A full 
three-dimensional numerical solution of a closely related 
situation was recently given || in connection with exper- 
imental work at NIST ||. Both of these analyses rely 
on the Gross-Pitaevskii equation, which effectively as- 
sumes a Hartree ansatz with a fixed number of particles, 
or spontaneous symmetry breaking. In addition, a recent 
analysis of matter-wave phase conjugation treats the cen- 
tral mode to all orders, but describes the side-modes in 
a linearized fashion j^]. In contrast, the present theory, 
which also holds for small condensates, does not make 
any of these assumptions. It handles all modes on equal 
footing, to all orders, and makes no assumption about 
their statistical properties. We note that our approach 
is closely related in spirit to a recent analysis of matter- 
wave four- wave mixing Q| , which is also valid to all orders 
and uses angular momentum algebra to diagonalize the 
problem. 
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II. PHYSICAL MODEL 



V f = ^2\f,m)(f,m\ 



(8) 



We consider a Bose-Einstein condensate of 23 Na atoms 
in the F = 1 hyperfine ground state, with the three in- 
ternal atomic states \F = l,mp = — 1), \F = 1, mp = 0) 
and |F = l,m.p = 1) of degenerate energies in the ab- 
sence of external magnetic fields. The condensate is con- 
fined by a far-off-resonance optical dipole trap. It is de- 
scribed by the three-component vector Schrodinger field 

*(r,t) = {*_ 1 (r,t),^ ( r ,i),*i(r,i)} (1) 

which satisfies the bosonic commutation relations 

[* i (r,t),tf j (r',t)]=6 ij 6(r-r'). (2) 

Accounting for the possibility of two-body collisions, its 
dynamics is described by the second-quantized Hamilto- 
nian 

H = J rfr*t(r,t)7J *(r,i) 

+~ J dr 1 dr 2 *t( ri)t )^t( r2)f )y( ri _ r2 )^( r2)t )^( ri)t ) ) 

where the single-particle Hamiltonian is given by 

H = p 2 /2M + V trap (4) 

and the trap potential is of the general form 
+i 

V trap = ]T U{v)\F=\,m){F = l,m\. (5) 

171= — 1 

Here p is the center-of-mass momentum of the atoms of 
mass M and £/(r), the effective dipole trap potential for 
atoms in the 1 1 , m) hyperfine state, is independent of m 
for a non-magnetic trap. 

The general form of the two-body interaction V(ri — 
r 2 ) has been discussed in detail in Refs. We repro- 

duce its main features for the sake of clarity. Considering 
exclusively situations where the hyperfine spin Fi = 1 of 
the individual atoms is preserved, we label the hyperfine 
states of the combined system of two collision partners 
with total hyperfine spin F = Fi + F2 by \f,m), where 
/ = 0, 1, 2 and m = — /, ...,/. In the shapeless approxi- 
mation, it can be shown that the two-body interaction is 
of the general form M 



where 



and 



V(r t - r 2 ) = S(ri - r 2 ) ^ hgfVf, 



.<?/ = Airhcif /M, 



(G) 



(7) 



is the operator which projects the state of the atomic 
pair onto a state of total hyperfine quantum number /. 
Here a/ is the s-wave scattering length for the channel 
of total hyperfine spin /. As a result of the symmetry 
requirement for bosonic atoms, it can be shown that only 
states with even / contribute to V(ri — r 2 ), so that 

V(r 1 -r 2 ) = M(ri - r 2 )(g 2 V2 + .9oPo) 

= fttf(n - r a ) (co + c 2 Fr • F 2 ) . (9) 



In this expression, 



co 
c 2 



(50 
(52 



2 ff2 )/3 
#o)/3, 



(10) 



which follows from the identities V\ + V2 = I and 
Fi • F 2 = V 2 - 2"P - Substituting this form of the 
two-body potential y(i"i — r2) into the second-quantized 
Hamiltonian (0) leads to 



(3) ^_ 



n = J2j ^L(r,t) 



— + U(r) 
2M K 1 



*™(r,<) 



+- / dr{{c +c 2 )[* t 1 * t 1 * 1 * 1 + * t _ 1 * t _ 1 *_ 1 *_ 1 

+2*J* (*l*i + *Li*-i)] + c *S* *o*o 
+2(c - ca)*!*!* 1 ^*-! 

+202(^*^*0*0 + ^.)}. (11) 

The physical interpretation of the various terms of this 
Hamiltonian has been discussed previously: The three 
terms quartic in one of the field operators, i.e. those 
of the form *^*j*i*i are self-(de)focussing terms, the 
terms involving two hyperfine states conserve the popu- 
lations of the individual spin states and merely lead to 
phase shifts, and the terms involving the central mode 
*o and both side-modes correspond to spin-exchange col- 
lisions. This "four-wave mixing" interaction, involving 
e.g. the annihilation of a pair of atoms with rap = and 
the creation of two atoms in the states rap — ±1, leads to 
phase conjugation in quantum optics and to matter-wave 
phase conjugation in the present case S. 

The analogy with optical four-wave mixing becomes 
even more apparent when we consider a situation where 
atoms in the rap = state are placed in a linear super- 
position of two counterpropagating center-of-mass modes 
of momenta ±hko, that is 



* o (a0 = -4= (e lkax a 01 + e- ik ° x a 02 ) , 
V V 



(12) 



while the atoms of spin mp = ±1 are taken to be at rest 



(13) 
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Here aoi and 002 are the annihilation operators of the two 
counterpropagating mj? = modes, with [aoij%j] = 
i,j = 1 or 2, and ai,a_i are the corresponding opera- 
tors for the modes associated with mp = ±1. Finally, 
V is the is the confinement volume of the condensate. 
Inserting this mode expansion into the Hamiltonian ( |ll| ) 
and ignoring all non-phase matched contributions finally 
yields the four-wave mixing Hamiltonian 

h 2 kn - hen ~ , - 

H = 2-iV + —1n(N - 1 

2M 2 V ; 

hC2 , + + ++ r< t t 

H — (a{a{aiai + a_ 1 al_ 1 a— ia_x — ialaia^a-i 



[JV",W] = 



+ 2(a\a 1 + a t _ 1 a^ 1 )(al 1 a al + aj 2 a 2) 
+ 4a}at 1 aoiao2 + 4aia_iaJ 1 a ( ' )2 ), 



(14) 



where we have introduced the total number of atoms 

N = a\ai + a^jO-i + aJjOcii + aj 2 a 2. (15) 

Note that the coupling between the central and side- 
modes of the condensate requires both energy and mo- 
mentum conservation. Energy conservation can be 
achieved e.g. by introducing an external dc magnetic 
field. Then the quadratic Zeeman effect compensates the 
mismatch in kinetic energies and provides the necessary 
condition E\ + E-\ = % k^/m (E±l being the energies 
of mp — ±1 spin states ). Indeed, for the magnetic fields 
used in Ref. |) E x + E- X ~ 1(T 28 J which is of the same 
order as the central mode kinetic energy h 2 k 2 /2m for 
fc ~ 10 7 m. 



In addition, the population 



erences 



Do — OQl a 01 — Qo2 a 02 

D = a\ai — a 2 a2 



(16) 



(17) 
(18) 



are also conserved. These two conservation laws follow 
directly from the fact that the Hamiltonian ( [hi]) describes 
the creation and annihilation of bosonic atoms in pairs. 
For example, the annihilation of two atoms in the ccntcr- 
of-mass modes 1 and 2 of the mp = hyperfine state, 
described by the operator pair 001002, results in the cre- 
ation of a pair of atoms in the mp = ±1 spin states via 
a\a i _ 1 . 

Together with the conservation of the total number of 
atoms, these conservation laws yield the two additional 
conserved quantities 



and 



= a\a\ + a 01 a i 



N 2 = a!]a_i + aJ 2 ao2- 



(19) 



These conservation laws make it possible to define an 
angular momentum algebra for this four-wave mixing sys- 
tem analogous to the Schwinger coupled boson represen- 
tation used in Ref. |llj for the description of two-mode 
condensates and in Ref. 0] for the three-mode coupling 
problem. However the present situation requires that 
these considerations be extended to a compound angular 
momentum representation Jl^ ]. 

We proceed by introducing spinor operators ai and &2 , 



III. ANGULAR MOMENTUM 
REPRESENTATION 



ai 



aoi 
«i 



and a2 



o 2 

O02 



(20) 



Despite the formal analogy between the Hamiltonians 
describing optical and atomic four- wave mixing, there are 
important details in which these processes differ. One of 
them is phase matching, which is normally more difficult 
to achieve with de Broglie waves due to the quadratic 
dispersion relation of atoms. However, the use of multi- 
component condensates permits to avoid this difficulty, 
as we have just seen. More important perhaps is the fact 
that nonlinear optics experiments usually use coherent 
laser light as pumps. It is well known that the num- 
ber of photons in these fields is not well determined. In 
contrast, the number of atoms in a condensate is a fixed 
(and integer) quantity in the absence of losses. This leads 
to important conservation laws, that we exploit in the 
following to reexpress the Hamiltonian (l3) in terms of 
angular momentum operators. This leads in turn to an 
exact solution of the problem. 

The most obvious conserved quantity of the Hamilto- 
nian ([l5]) is the total number of atoms, as follows from 
the commutator 



as well as the two angular momentum operators Si and 
S a , 



I ( a 01 a l + OOlOl 
Si = ajerai = - -i{a\ x ax - Ooiaf) 
1 \ t t 



(21) 



and 



S 2 = a^cra2 



a 2°02 + a Q2 a 2 



^ 1 "-2"uz ' "-02"-* \ 

- -i(a 02 al - al 2 a 2 ) , (22) 



a\a,2 — oj 2 ao2 



where <r is the Pauli spin operator. The Casimir opera- 
tors Kj associated with the angular momenta Sj are also 
constants of motion since 



Expressed in terms of the total spin operator 



(23) 
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Si +s 2 



(24) 



with components Sj 
becomes 



U 



2M 



2hc? 



- Sij + S 2 j, the Hamiltonian ([jl 




(25) 



Together with the conserved Casimir operators, this form 
of the four-wave mixing Hamiltonian suggests expressing 
the state of the system in terms of the complete set of 
quantum numbers associated with the eigenstates of the 
operators S 2 , S 2 , S 2 and S z . The Hamiltonian ( p5|) can 
clearly be diagonalized in terms of these states, yielding 
a complete set of eigenstates and eigenenergies. 

To illustrate how this works, we consider for concrete- 
ness a condensate consisting of N atoms such that iti- 
tially iVi = N 2 = N/2, to < N/2 atoms in the hyperfine 
state m,F = 1 and none in the state mf = 0. In the 
"natural" eigenbasis of the operators {K\, S\ z , K 2 , S 2z } 
the initial state of this system is described by the state 
vector 



10(0)) = 



N N 



N 



4 1 



(26) 



Expanding it in terms of a complete set of eigenvectors 
{|f , f , S, S z )} of the Hamiltonian @, we find 

s 

(27) 



where 
C(S lz ,S 2z ;S, -to) 



,AT TV 
M 'T 



-m — , Si z ; —,b 2z ) 



(28) 



are Clebsch-Gordan coefficients |L3| |, which are nonzero 
only for S z = S\ z + S 2z = -m. 

This conservation of S z under the Hamiltonian fl25| ) 
further implies that the state vector of the system is given 
at time t by 



Wt))=X> s (t)l T>T ;£- 



s 



as(t)\S, -to) 



(29) 



where in the second equality we have made explicit use 
of the value of the conserved quantities K\ , K 2 and S z 
to simplify the notation via 



N N 

1 



S, —to) — > jS 1 , —to). 



(30) 



Note that this simplification is not general, but is appro- 
priate for the initial condition at hand. 

The equations of motion for the probability amplitudes 
cts(t) follow from the Schrodinger equation. For the con- 
densate containing N atoms they read, in a frame ro- 
tating at the frequency hk 2 N/2M + c N(N ~ l)/2 - 
2c 2 (N/2 + (D /2) 2 ), 



ia s (jk) = 2c 2 S{S + l)as(t) 



(31) 



and can be integrated trivially with the initial condition 

N N 
as{t = 0) = C(— - to, - — ; 5, -to). 



IV. DYNAMICS 

First we apply the results of the preceding section to 
the study of the dynamics of population exchange be- 
tween the different modes of the condensate. For in- 
stance, the population of the mf = 1 hyperfine spin state 
can readily be determined from the expectation value of 
Siz- With Eqs. ([l9|) and the definition of Sx z we find 



(a\ai) 



N 



(Si 



From Eq. (|29j) we have 

{Siz) = ^ot s {t)C(p,-{p + m);S,- 

P s 



where we inserted the identity operator 
I 



E 

Pl,P2 



N N ,,N N 
T,Pi; T,P2)(-r,Pi; -riVl\ 



(32) 



(33) 



(34) 



and used the simplified notation Q3y) as well as the prop- 
erty that the Clebsch-Gordan coefficients' are nonzero 
only for S z — S\ z + S 2z . Eq. ([53]) can be evaluated 
numerically. 

The evolution of the population of the mp = 1 side- 
mode is shown in Fig. 1 for N — 100 atoms in the system. 
In case (b) the initial mode population is (a\a\) = m = 5, 
while case (a) illustrates the build-up from noise, m = 0. 
In both cases, the sidemode population exhibits an ini- 
tial growth to the point where it contains about 1/3 of 
the atoms in the first case and about half of the atoms 
in the second. This is followed by a collapse to a quasi- 
steady state population, as well as a subsequent revival 
at 2c 2 t\ = 7r. This dynamics then repeats itself period- 
ically, with revivals at 2c 2 t n = nn, independently of N. 
This is similar to the periodical revivals which occur in 
two-photon Janes-Cummings model discussed in fl4|| . 

During the periods of collapse, one has Sj Z = —to/2, so 
that all modes are almost equally macroscopically popu- 
lated with (aJ x aoi) = N/4 - to/2, {a\ai) = N/4 + to/2, 
( a J 2 a 02 ) = N/4 + m/2 and {a\a 2 ) = N/4 - to/2. 
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A particularly interesting aspect of the present study 
is that it allows one to obtain the quantum correlations 
between sidemodes. In optics, for example, four-wave 
mixing provides a method to study purely quantum me- 
chanical effects such as squeezing and nonclassical states 
of the radiation field, and also to prepare states of com- 
posite systems exhibiting strong quantum mechanical en- 
tanglement |l5j . These states are of considerable interest 
in tests of the foundations of physics as well as quan- 
tum information processing such as quantum cryptog- 
raphy [|l6],[l7| and quantum computing |L8j . Macroscopic 
quantum states of massive particles present an interesting 
alternative to all-optical systems, hence it is of consider- 
able interest to determine to which extent quantum en- 
tanglement between sidemodes can be achieved in Bose- 
Einstcin condensation. 

In analogy with the optical case, one can quantify the 
amount of quantum entanglement between condensate 
modes by determining the extent to which the Cauchy- 
Schwartz inequality is violated by the second-order cross- 
correlation functions between modes fl5|] . In particular, 
for a 'classically looking ' optical system with positive 
Glauber P-representation, the single-time second-order 
cross-correlation function is bound by 



where 



G%(t)<[G\ 2) (t)Gf\t)}V 2 . 



(35) 



In case the "P-representation is not positive or does not 
exist, in contrast, the upper bound is higher, namely 



Gg (t) < { [Gf (t) + Gf (i)) (Of (t) + Gf (t)] } 



1/2 

(36) 



In these inequalities, we have introduced the single-time 
and single-mode first-order correlation functions 



Gf\t) = m)\ahm)) 



(37) 



as well as the single-time two-mode second-order corre- 
lation functions 

G§\t) = {m4^hMt)) (38) 
and the single-time, second-order correlation functions 

Gf\t) = {m4 a h a i\<t>it))- (39) 

As was the case for the sidemode populations, the 
single-time single-mode second-order cross-correlation 
between the sidemodes uif = ±1 can be expressed in 
terms of the z-component of the individual pseudospins 
as 



<o\ N N 

G ( %(t) = ((--S lz )(--S 2z )) 

= (t) ~ T {{Slz) + {S2z)) + {SlzS2z) > (40) 



= ^2p(-P- m) 



[N/2] 

^2 a S (t)C (p, -{p + m);S, 



m , 



(41) 



and the sum can easily be evaluated numerically. 

Figure 2 compares the time dependence of the normal- 
ized central-mode — side-mode correlation function (lower 
curve) 



(42) 



^G&\t)G[ 2 \t) 
and the side-mode — side-mode correlation (upper curve) 



G%{t) 



G {2l (t)G^{(t) 



(2), 



(43) 



for (a) the case m = where the sidemode builds up 
from quantum fluctuations and (b) the case m = 5 of an 
injected signal. These results illustrate how the corre- 
lations between central-mode and side-modes do satisfy 
the classical Cauchy- Schwartz inequality while the side- 
mode — side-mode cross correlations violate them. The 
violation is particularly strong in the case of build-up 
from noise, as should be intuitively expected. In that 
case the hyperfine sidemodes rap = ±1 play symmetric 
roles, thus 



G?L(t) = GHt) + GT(t) 



(44) 



i.e. Eq. ( pq ) thus becomes an equality, corresponding to 
the maximum violation of the classical Cauchy-Schwartz 
inequality allowed by quantum mechanics. 

The difference in the behavior of the two-mode corre- 
lation functions between side-modes and those involv- 
ing one central and one side-mode can be intuitively 
understood from the form of the wave-mixing term 
a\a ! _ 1 aoiao2 appearing in the Hamiltonian (|l4|). Indeed, 
the coupling between side-modes, involving two annihi- 
lation operators, is reminiscent of the interaction a\a 2 in 
the Hamiltonian of parametric amplification leading to 
squeezing and quantum entanglement between two side- 
modes. In contrast, the coupling between central and 
sidemodes involves both an annihilation and a creation 
operator. 



V. SUMMARY AND CONCLUSIONS 

In contrast to the optical case, where it is usually diffi- 
cult to study quantum mechanically four-wave mixing in 
the strong field limit, this task can be achieved relatively 
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easily for matter waves, a direct consequence of the con- 
servation of the total number of atoms in a condensate at 
T = and in a lossless trap. The resulting conservation 
laws permit to develop an angular momentum algebra 
analysis that leads to an exact solution of the problem 
away from the linear regime where the sidemode popula- 
tions remain small. 

In this paper we have applied this technique to the 
study of the dynamics of the population exchange be- 
tween hyperfine levels of a condensate. We found that 
this exchange is periodic and it is characterized by a se- 
quence of "collapses" and "revivals" reminiscent of those 
appearing in the two-photon Jaynes-Cummings model 
of optical physics. In addition, we found that strong 
quantum correlations can develop between the central 
modes and the side modes, the general state of the sys- 
tem exhibiting a strong quantum entanglement between 
the modes rap = ±1. Thus, it appears that multicom- 
ponent condensates offer a fascinating method to create 
quantum entanglement at a truly macroscopic level, a 
possibility made even more attractive by the fact that 
these systems suffer very little from dissipation, since 
they consist of ground-state atoms. The finite lifetime 
of condensates is usually attributed to three-body colli- 
sions, which result in losses on a timescale of seconds. 

We conclude by noting that the revivals in the popula- 
tion exchange occur at a time independent of the number 
of atoms in the condensate. Hence, they allow for a di- 
rect and absolute determination of the coefficient C2. In 
practice, however, C2 seems to be too small to observe 
revivals on the timescales of the condensate lifetime. 
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FIG. 2. Time evolution of the one-time normalized 
central-mode — side-mode correlation function (lower curve) 
TZ^i 1 (t) and the one-time normalized side- mode — side-mode 
correlation function (upper curve) TZ^\ t (t) for N=100 and 
(a)m=0, (b) m=5 
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